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Recently, stereotactic body radiotherapy (SBRT) for lung cancer is conducted with heterogeneity-corrected 
treatment plans, as the correction greatly affects the dose delivery to the lung tumor. In this study, the correl- 
ation between the planning target volume (PTV) and the dose delivery is investigated by separation of the het- 
erogeneity correction effects into photon attenuation and electron transport. Under Institutional Review Board 
exemption status, 74 patients with lung cancer who were treated with SBRT were retrospectively evaluated. 
All treatment plans were generated using an anisotropic analytical algorithm (AAA) of an Eclipse (Varian 
Medical Systems, Palo Alto, CA) treatment planning system. Two additional plans were created using the 
same treatment parameters (monitor units, beam angles and energy): a plan with no heterogeneity correction 
(NC), and a plan calculated with a pencil beam convolution algorithm (PBC). Compared with NC, AAA and 
PBC isocenter doses were on average 13.4% and 21.8% higher, respectively. The differences in the isocenter 
dose and the dose coverage for 95% of the PTV (Dgs%) between PBC and AAA were correlated logarithmically 
(p = 0.78 and p = 0.46, respectively) with PTV. Although 095% calculated with AAA was in general 
2.9% larger than that for NC, patients with a small PTV showed a negative ADgs% for AAA due to the signifi- 
cant effect of electron transport. The PTV volume shows logarithmic correlation with the effects of the lateral 
electron transport. These findings indicate that the dosimetric metrics and prescription, especially in clinical 
trials, should be clearly evaluated in the context of target volume characteristics and with proper heterogeneity 
correction. 
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INTRODUCTION 

Stereotactic body radiotherapy (SBRT) has gained world- 
wide acceptance as treatment for early stage inoperable small 
lesions of non-small cell lung cancers (NSCLCs), and has 
yielded excellent tumor control rates, above 90% [1, 2]. 
Several clinical trials have been conducted to investigate the 
use of SBRT for NSCLC treatment. Radiation Therapy 
Oncology Group (RTOG) protocols 0236 and 0618 did not 
allow heterogeneity corrections for dosimetry and monitor 
unit calculation because of a lack of clinical experience in 
heterogeneity correction. Later, RTOG protocols 0813 and 
0915 did require heterogeneity corrections in dose calcula- 
tion [3]. Dose calculation in lung SBRT is complicated as a 



consequence of the dosimetry of thoracic organs with tissue 
heterogeneities and small treatment fields [4]. The accuracy 
of dosimetry has improved a great deal, with advanced com- 
mercial dose calculation algorithms such as collapsed-cone, 
convolution/superposition, anisotropic analytical algorithm 
(AAA), Acuros® XB and Monte Carlo [5-9]. The Japan 
Clinical Oncology Group (JCOG) 0403 phase II trial and 
0702 phase I trial were conducted using heterogeneity cor- 
rection with one-dimensional equivalent path length (EPL) 
and a convolution/superposition algorithm, respectively [10]. 
Recently, some clinical trials have been conducted using 
advanced dose calculation algorithms. 

Radiation dosimetry plays an important role in the com- 
parison of clinical outcome. In light of this, the impact of 
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heterogeneity corrections on dose calculation in lung SBRT 
has been investigated [11-14]. Many of the earlier studies 
that investigated the effect of heterogeneity corrections using 
clinical data assessed a small number of cases or collected 
patient data from multiple institutions. Xiao et al. [15] and 
Schuring and Hurkmans [16] reported that the dose covering 
95% of the PTV (£>95%) for a PTV calculated with superpos- 
ition or collapsed-cone was smaller than that for a PTV cal- 
culated without heterogeneity corrections. Recently, Ueki 
et al. [17] investigated 83 SBRT plans and reported the op- 
posite for Z)95% results. In their study, the D()s% of a PTV cal- 
culated with AAA was slightly larger than that calculated 
without corrections, however they did not mention the dis- 
crepancy between previous reports and their study. The 
effect of heterogeneity corrections on target coverage is crit- 
ical when determining treatment protocols. Although previ- 
ous studies evaluated the average ± SD of the dose-volume 
metrics, or statistical differences among various dose calcula- 
tion algorithms, few studies have investigated the factors that 
affect the heterogeneity corrections. 

To evaluate more detailed effects of heterogeneity correc- 
tions with advanced algorithms on the target dose of SBRT 
for lung cancer, a retrospective analysis of the treatment 
plans of 74 patients was conducted at a single institution. 
The results were analyzed according to size of target volume 
and by separating photon attenuation and electron transport 
for the heterogeneity corrections. 

MATERIALS AND METHODS 

Patients 

Under Institutional Review Board (IRB) exemption, we 
retrospectively analyzed the treatment plans of 74 patients 
with lung cancer who were treated with SBRT at Indiana 
University School of Medicine, Indianapolis. There were 22 
and 52 cases in the right- and left-lobe of the lung, respect- 
ively. Details of the prescribed doses for these patients are 
listed in Table 1. Every SBRT patient was immobilized 
using either an Elekta Stereotactic Body Frame (Elekta, 
Stockholm, Sweden) or a CIVCO Body Pro-Lok system 
(CIVCO, Kalona, lA), depending upon the patient size, 
comfort and suitability as decided at the time of simulation. 

Treatment planning and evaluations 

Eclipse version 10.0 was used for the treatment planning. 
Plans were generated to cover 95% of the PTV with the pre- 
scribed dose, thus providing typically a 25% higher dose to 
the GTV. Eclipse provides two algorithms for heterogeneity 
correction; pencil beam convolution (PBC) and AAA. The 
PBC model in Eclipse is mainly for homogenous medium 
(water) dose calculation that gets supplemented with older 
models modified Batho, Batho power law and equivalent 
tissue air ration (ETAR) [18]. It was noted that both modified 
Batho and Batho power law gave almost the same results for 
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Table 1. Description of treatment parameters 

Plan parameters Frequency 

Energy 

6 MV 50 (68%) 

e+ieMV" 24(32%) 
Number of beams 

8 3 (4%) 

9 15 (20%) 

10 49 (66%) 

11 7(10%) 
Number of non-coplanar beams 

<2 2 (3%) 

4 18 (24%) 

5 46 (62%) 

6 8(11%) 

"6+16 MV = treatment fields including both 6 MV and 16 MV 
beams. 

SBRT, and hence default modified Batho with PBC was 
used. The ETAR option with PBC cannot be used with non- 
coplanar field arrangements, thus ETAR was not an option in 
this study as the majority of the SBRT fields were non- 
coplanar. Heterogeneity correction was applied to all clinic- 
ally approved plans using an AAA algorithm. The AAA is 
considered a superior algorithm compared with older and 
pencil beam algorithms [7, 19-21]. Two additional plans 
were generated for each patient: (i) a treatment plan with no 
coiTection (NC), and (ii) PBC with modified Batho hetero- 
geneity correction. In the PBC algorithm, the dose deposited 
at a point was calculated as a convolution of energy fluence, 
or total energy released per unit mass (TERMA), with the 
respective dose deposition kernel pre-calculated for a narrow 
beam in water [22]. The Batho power-law correction method 
is an empirical correction to account for both primary beam 
attenuation and scatter changes in heterogeneous materials. 
The modified Batho correction uses only the descending part 
of the TAR/TMR curve because the curve in the build-up 
region of a high-energy photon is no longer valid. However, 
PBC does not take into account changes in lateral electron 
transport. This provided us with an opportunity to differenti- 
ate between photon attenuation (PBC) versus electron 
transport (AAA). 

For each patient, two additional plans (NC and PBC) were 
created using the same monitor units as calculated with 
AAA. Because PBC considers changes in tissue density but 
does not consider lateral electron transport, the difference 
between the NC and PBC can be considered as the effect of 
photon attenuation. As AAA considers both depth correc- 
tion and electron transport, the difference between PBC and 
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AAA therefore represents the effect of electron transport 
calculation. 

To evaluate the effect of photon attenuation, we assessed 
the change in path length (APL) defined as: 



(1) 



where dj and EPLj represent the physical and equivalent path 
length (EPL) depth for each treatment beam, respectively. Wj 
and W represent the weight of each beam and the total 
weight, respectively. The Dgs% and the volume receiving the 
prescribed dose (yioo%), representing the quality of the target 
volume coverage, were also analyzed. 

Statistical analysis 

JMP software (ver. 9.0.2; SAS Institute, Gary, NC) was used 
for statistical analysis. All pairwise comparisons among the 
three calculation algorithms were conducted for the isocenter 
dose (Z)iso)> ^95% and Vioo% using the Steel-Dwass test. The 
correlation between dosimetric parameters and anatomical 
characteristics including APL, the PTV and the distance 
between the PTV and the chest wall or mediastinum was 
assessed with Spearman's rank correlation coefficient. 
Statistical significance was defined as a /'-value < 0.05. 

RESULTS 

The characteristics of target volumes of the patients are 
shown in Table 2. Due to the selection criterion in our insti- 
tution, peripheral lesions were preferred for SBRT, as shown 
in Fig. 1, which indicates the frequency distribution of 
lesions. This figure shows that the majority of the patients 
(74.3%) have the GTV close to the chest wall (<0.5 cm) and 
the remainder have it distributed evenly up to 2.5 cm. 
Figure 2 shows an example of dose distributions using the 
three treatment plans (NC, PBC, AAA). The subtle dose 
difference is clearly visible. Typically, PBC generally overes- 
timates dose in lung as it underestimates the range of the sec- 
ondary particles, whereas AAA provides accurate dose 
distribution. 



Table 2. Details of target volume 



Anatomic parameters 



Mean ± SD (min - max) 



Target volume 




GTV 


13.8 ±16.5 (0.6-1 11.5) [cm^] 


PTV 


45.3 ±35.7 (9.5-212.0) [cm^] 


Maximum dimension 




GTV 


3.34 ± 1.31 (1.26-7.02) [cm] 


PTV 


5.12± 1.33 (1.54-8.75) [cm] 



Isocenter dose 

The Diso relative to the prescribed dose for each plan is 
shown in Fig. 3a. The Djso (average ± SD) for the entire 
patient population was 110.8 ±4.4%, 132.6 ±4.3% and 
124.2 ±3.1% for NC, PBC and AAA, respectively. The 
Z^iso calculated with PBC and AAA were larger than that 
calculated with NC by 21.8% and 13.4%, respectively, in 
our population, with both P < 0.0001. The Diso of AAA was 
8.4% smaller than that of PBC (P < 0.0001). In Fig. 3b, the 
isocenter dose difference (ADjso) between the dose with and 
without heterogeneity corrections of PBC (APBC) and AAA 
(AAAA) were plotted with APL as shown in Eq 1 . The slope 
of the lines for PBC and AAA are nearly identical. The 
APBC and AAAA showed intermediate linear correlation 
with APL (p = 0.60, P < 0.0001 for both). 

Target coverage 

Typically, patients are treated with better than 95% target 
coverage, as shown in Fig. 4a indicating Vioo% with NC, 
PBC and AAA algorithms for all patients. The median Vioo% 
was 90.6% (range, 42.3-98.7%), 100.0% (range, 92.9- 
100.0%) and 96.0% (range, 75.6-99.9%) for NC, PBC and 
AAA, respectively. The long error bars indicate the variability 
among the patients. The average ± SD of 1)95% was 
97.6 ±5.0%, 114.6 ±7.4% and 100.5 ±2.8% for NC, PBC 
and AAA, respectively as shown in Fig. 4b. The 1)95% of 
PBC and AAA were larger than that of NC by 17.0% 
(/■< 0.0001) and 2.9% (/'< 0.0001), respectively. The 
ADi)s% of AAA plans was 14.1% smaller than that of PBC 
(P < 0.0001). Figure 4c illustrates the evaluation of ADgs% of 
the PTV for each patient, representing a similar analysis to 
that of Fig. 3b for Diso- The APBC shows intermediate linear 
coiTelation with APL (p = 0.58, P < 0.0001), although AAAA 
showed weaker correlation with APL (p = 0.41, P = 0.0003). 

Correlation with target volume 

The size of the PTV is correlated with ADiso and ADgs% 
between the dose with and without heterogeneity corrections. 
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Fig. 1. Frequency distribution of the shortest distance (8) between 
gross tumor volume (GTV) and chest wall or mediastinum. 
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Fig. 3. (a) The isocenter dose (Disa) calculated without (NC) and 
with (PBC with modified Batho and AAA) heterogeneity 
corrections are shown. Boxes represent median and quartile values. 
Whiskers represent maximum and minimum values, (b) The 
differences between Dis„ calculated without and with heterogeneity 
corrections (ADiso) are illustrated. 



Negative distance represents the largest overlapped distance 
between the PTV and the chest wall. Although isolated 
tumors with positive distance showed moderate linear correl- 
ation with the distance (r^ = 0.62, /'< 0.0001), whole cases, 
including tumors attached to chest wall, showed weaker 
coiTelation (r^ = 0.34, f'<0.0001). 

In Fig. 6, the differences in the dose distribution have been 
illustrated for representative patients with (a) large (102.4 
cm^) and (b) small (13.2 cm^) PTV volumes. The APL were 
3.5 cm and 0.8 cm for patients (a) and (b), respectively. 
Patient (a) showed a larger attenuation effect in the APBC as 
a consequence of a large lung:solid tissue ratio in the path 
length. In the A(AAA - PBC) of patient (b), the larger lateral 
electron transport effect is visualized, resulting in the nega- 
tive net effect shown in AAAA. 

DISCUSSION 



as shown in Fig. 5a and b, respectively. The dose differences 
between PBC and AAA A(AAA-PBC), are also plotted. 
Interestingly, the ADjso of A(AAA-PBC) showed sig- 
nificant logarithmic relationship with the PTV volume 
(p = 0.78, P < 0.0001). Almost all AAAA cases showed posi- 
tive ADiso, although one patient showed a negative value. 
This result indicates that the effect of attenuation was larger 
than that of electron transport for most cases. As illustrated 
in Fig. 5b, the ADg5% of A(AAA - PBC) showed intermedi- 
ate logarithmic correlation with the PTV volume (p = 0.46, 
/"< 0.0001), although the result showed a larger deviation 
than that of ADiso- Most cases of ziAAA showed a posi- 
tive value of ADg^%, but for a small PTV volume it showed 
negative values, probably as a consequence of significant 
effects of electron transport. In Fig. 5c, the ADgs% of 
A(AAA-PBC) were plotted against the shortest distance 
between the PTV and the chest wall or mediastinum. 



Several studies have investigated the impact of heterogeneity 
corrections on the dose distribution in lung SBRT with 
similar findings. We demonstrated that the £>95% of a PTV 
calculated with AAA is 2.9% larger than that of a PTV calcu- 
lated with NC (Fig. 4b). Interestingly, this is the opposite of 
£)g5% results reported by Xiao et al. (-1.6%) and Schuring 
and Hurkmans (-5.2%), although our results of Di^o 
(13.4%) were similar to the previous report by Xiao et al. 
(12.5%) [15, 16]. Recently Ueki et al. investigated 83 
patients and reported that the PTV 1)95% under AAA calcula- 
tion was higher than that without corrections by -1.4% [17]. 
Their results for 0^^% were similar to ours, although they 
did not investigate the factors that affect the heterogeneity 
coiTections. 

We demonstrated the significant logarithmic relationship 
between the effect of lateral electron transport and PTV 
volume (Fig. 5a and b). Cases with a large PTV showed 
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smaller effects from the lateral electron transport than those 
from attenuation, resulting in the positive ADgs% of AAAA. 
In contrast, with a small PTV, the effects of electron transport 
become larger (Fig. 6b). Although the average PTV 0^5% of 
AAA plans is 2.9% larger than that of NC, it depends on the 
anatomical characteristics of the patients analyzed in this 
study. Schuring and Hurkmans previously investigated the 
effects of heterogeneity corrections on the dose distribution 
of 26 SBRT patients and reported that the collapsed-cone 
convolution algorithm showed decreased conformity with 
the PTV volume [16]. On the other hand, van der Voort van 
Zyp et al. compared EPL and the Monte Carlo algorithm 



with respect to tumor location and size [23]. They separated 
patients roughly by tumor size and reported that the differ- 
ence in Dg^% between two algorithms was larger for small 
tumors. The results of previous studies that investigated 
various calculation algorithms cannot be simply compared 
without considering the variation in target volume character- 
istics. The lateral electron transport is affected by the 
medium density, size of the target volume and also by the 
surrounding tissue, such as chest wall. Narabayashi et al. 
[24] investigated the lung SBRT plans for 20 patients and 
reported that the ratio of MU values calculated with the 
Batho Power Law and Monte Carlo showed linear correlation 
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Fig. 6. Differences in dose distribution for APBC, AAAA and A(AAA - PRC). The representative patients with (a) large, and 
(b) small tumors are illustrated. 



with the distance between the PTV and the chest wall. We 
also analyzed the correlation between the effects of lateral 
electron transport and the distance between the PTV and the 
chest wall or mediastinum (Fig. 5c). Although the data of 
isolated tumors showed moderate linear coiTelation, the 
whole cases (including tumors attached to the chest wall) 
showed weaker correlation. To properly assess the effect of 
the surrounding tissue on the lateral electron transport, the 
area of attachment to the surrounding tissue, and the volume 
of the chest wall included in the PTV (and its density) should 
also be considered. These complexities of the lateral elec- 
tron transport effects caused by the surrounding tissue will 
lead to larger Dgs% variation of A (AAA - PBC), illustrated 
in Fig. 5b. 

In the current study, the treatment plans were clinically 
optimized based on AAA; and the beam arrangements, 
energy and MU were not modified when calculating with 
NC and PBC. The differences in dosimetric parameters 
between AAA and the other two algorithms become smaller 
with optimizing plans, although the results include the uncer- 
tainty due to the individual planning techniques. Some 
studies compared commercial Monte Carlo dose calculation 
with other algorithms for dose calculation of lung SBRT [3, 
25, 26]. Li et al. noted that some of the criteria of treatment 
protocols for supeiposition might be too strict for the Monte 
Carlo calculation to satisfy the requirements [3]. Several pre- 
viously reported that the advanced dose calculation algo- 
rithm showed worse target coverage. However, we showed in 
this study that the target coverage is strongly affected not 
only by calculation algorithms but also by the anatomic 
characteristics of the patients. 



We demonstrated the effects of the heterogeneity cor- 
rection by separating them into the attenuation and lateral 
electron transport components. The PBC provides photon at- 
tenuation, whereas AAA accounts for the electron transport 
that is needed for lung SBRT. The PTV volume shows sig- 
nificant correlation with the effects of the lateral electron 
transport with logarithmic coiTelation. When the patients are 
not classified with target volumes, the dosimetric parameters 
evaluated with average ± SD may not be able to assess the 
effects of the heterogeneity corrections properly. Several 
clinical trials for lung SBRT are being conducted using the 
target coverage-based prescription with heterogeneity correc- 
tions. It is concluded, therefore, that when radiation outcome 
is compared between clinical trials, dosimetric metrics, the 
target volume characteristics and heterogeneity correction 
algorithms should be properly accounted. 
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